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Abstract 

A continuum dislocation model of formation of grains whose boundaries have 
a non-vanishing thickness is proposed. For a single crystal deforming in sim¬ 
ple shear the lamellar structure of grains with thin layers containing disloca¬ 
tions as the geometrically necessary boundaries turns out to be energetically 
preferable. The thickness and the energy of this type of grain boundary are 
computed as functions of the misorientation angle. 
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1. Introduction 

One of the main guiding principles in seeking an appropriate theory of for¬ 
mation of grains in metals and alloys during and after cold working processes 
producing severe plastic deformations has hrst been proposed by Hansen and 
Kuhlmann-Wilsdorf [6] in form of the so-called LEDS-hypothesis: the dis¬ 
location structures in the hnal state of deformation minimize the energy of 
crystals (see also PDEIIII]). The main reason why the formation of grains 
becomes energetically preferable at severe plastic deformations lies in the 
non-convexity of the energy of crystal in this range [U m HSl E] . Within 
the conventional crystal plasticity considered in HEIE] the minimization 
of such non-convex energy leads immediately to the inhnitely hne lamellar 
structure with grain boundaries as sharp interfaces. However, as mentioned 
by Kuhlmann-Wilsdorf and Hansen HD], typical grain boundaries, termed 
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geometrically necessary boundaries, have as a rule a non-vanishing thick¬ 
ness and may contain a large number of dislocations and thus contradict the 
conventional crystal plasticity. The question then arrises in this connection: 
what kind of continuum model may resolve this conflict? The present paper 
proposes a dislocation model of formation of grains within the continuum 
dislocation theory IUDS] which predicts the existence of such geometrically 
necessary boundaries. By including the energy of dislocation network con¬ 
taining the gradient of the plastic slip into the energy functional we regularize 
the non-convex energy minimization problem. We illustrate the application 
of the theory on the example of single crystal deforming in single slip under a 
simple shear. We show that the geometrically necessary boundaries, in which 
the transition from one grain to the next occurs smoothly, have a small but 
hnite thickness and contain a large number of dislocations. Although the 
resultant Burgers vector of dislocations in such grain boundary is non-zero, 
they do not produce long range stresses, and the lamellar structure of grains 
is in fact the low energy dislocations structure. We also compute the thick¬ 
ness of geometrically necessary boundaries and their energies as functions 
of the misorientation angles. Based on these results we estimate also the 
number of grains in terms of the specimen sizes. 

2. Continuum theory of formation of grains 

We consider for simplicity an initially dislocation-free single crystal having 
only one active slip system. In this case the kinematic quantities character¬ 
izing its observable deformations are the placement held y(x) and the plastic 
slip held /3(x). The incompatible plastic deformation is given by 

FP(x) = I -I- /3(x)s ® m, 

with the pair of constant and mutually orthogonal unit lattice vectors s 
and m denoting the slip direction and the normal to the slip planes. Using 
the multiplicative resolution of the total compatible deformation gradient 
F = dy/dx into the plastic and elastic part [ 12 ], we hnd the incompatible 
elastic deformation in the form 

1 dy 

F^ = F ■ F^-^ = ^ ■ (I - ds 0 m). 
ox 

The tensor of dislocation density measuring the incompatibility of F^ reads 

(see [121 [ 13 ]) 

T = -F^ X V = s ® (V/? X m). 
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If, in addition, all dislocation lines are straight lines parallel to the unit vector 
1, then the scalar dislocation density (or the number of excess dislocations 
per unit area perpendicular to 1) can be determined as 

IT-11 1 

p=L^ = -|(V/?xm)-l|, 

with b being the magnitude of Burgers vector. 

For crystals having as a rule small elastic strains we propose the free 
energy per unit volume of the undeformed conhguration in the form 

V-(E',p) = b(trE')= +;,tr(E' • E') + (i) 

2 2 

Here E® = ■ F® — I) corresponds to the elastic strain tensor, A and 

/i are the Lame constants, k a material constant, and ps can be interpreted 
as the saturated dislocation density. The hrst two terms in (|^ represent 
energy of crystal due to the macroscopic elastic deformation. The last term 
describes energy of the dislocation network for moderate dislocation densi¬ 
ties. Note that for small or extremely large dislocation densities close to the 
saturated density the logarithmic energy proposed in [2] is more appropriate. 
We deform this crystal occupying in the initial conhguration some region V 
of three-dimensional space by placing it in a displacement-controlled device 
such that, at the boundary clV, the conditions 

y(x) = F ■ X, /?(x) =0 for X G aV (2) 

are specihed, with F being a given overall deformation. If the deformation 
process is isothermal, no body force acts on this crystal, and the resistance 
to the dislocation motion can be neglected, then the following variational 
principle turns out to be valid: the true placement vector y(x) and the true 
plastic slip /3(x) in the final equilibrium state of deformation minimize the 
energy functional 

^[y(x),/?(x)] = [ w(F,/3,VP)dx (3) 

Jv 

among all continuously differentiable helds y(x) and /9(x) satisfying con¬ 
straints (|^, where w(F,/3,V/3) = '0(E®,p). We will see that, due to the 
non-convexity of the free energy density ([^ and the presence of V/3 in the 
energy functional via the energy of the dislocation network, the formation of 
grains with regular grain boundaries having a hnite thickness is energetically 
preferable. 
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3 . Energy minimizer in plane strain simple shear 

Consider now the special case of plane strain simple shear of the specimen 
in form of a cnboid of height H, width L, and depth D snch that = x^, 
while yi (x), y 2 (x) and (3 (x) depend only on two cartesian coordinates Xi and 
X 2 and satisfy at the side bonndary the conditions 

yi = 7 x 2 , 2/2 = X 2 , (3 = 0 , 

with 7 being the overall shear strain. We assnme that s^ = (cos sin <^9,0), 
= (—sin cos 99 , 0 ) and all dislocation lines are parallel to the xa-axis, 
so p = |V/3 ■ s|/6. If the deformations are nniform snch that 

F = F = I + 761 (g) 02, F^ = I + / 3 s ® m, 

with 7 and j3 being the constants, then the energy ([^ normalized by p| V| and 
minimized with respect to f3 tnrns ont to be non-convex for p G (—7r/2,0) 
as shown in Fig. (see 0). 



Figure 1: Condensed energy 6 ( 7 ) and the dimensionless shear stress e'( 7 ) for ip = —45°. 

In view of this non-convexity, we proposed in [8] the energy minimizing 
seqnence consisting of layers having the nniform states A and B according to 


7A = 0, 

13a = Q, 


4 


7s = -2 cot (p, 
(3b = 2 cot (p. 






such that the volume fraction of the layer B is given by s = —7/(2 cot 99 ). 
It has been shown in [ 8 ] that such candidates for the minimizer satisfy the 
equations of equilibrium in each layer as well as the outer boundary condi¬ 
tions except at the side boundaries Xi = 0, L of the specimen. Besides, the 
energy of such lamellar structure is equal to zero which is the minimal pos¬ 
sible value. However, if the boundaries between layers are sharp interfaces, 
these candidates do not belong to the set of admissible helds of our origi¬ 
nal variational problem ([^ and (|^ due to the jumps of F and (3 on those 
interfaces, so they fail to be the energy minimizers of ([^ and (|^. 
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Figure 2: Enhanced grain boundary model. 

To correct the behavior of those candidates for minimizers we assume 
now that the layers corresponding to the states A and B are separated by a 
thin layer of small thickness h in which the placement and plastic slip change 
smoothly from state A to state B (see Fig. [^. Since this boundary layer is 
thin {h is much smaller than sif, (1 — s)H, and the sizes of the specimen), it 
is reasonable to assume that the displacement in X 2 -direction is zero, while 
the displacement in xi-direction and the plastic slip depend only on X 2 '- 


yi = xi + u{x2), y2 = X2, (3 = (3 {x2). 


With this Ansatz it is easy to show that the determination of functions u{x 2 ) 
and /3{x2) as well as the unknown boundary layer reduces to minimizing the 
following functional 

h[u{x2),/3{x2)]= [ [f {u^ 2,/3) + g{(3^2)]dx2 (4) 

Jh 
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among functions u{x2), ( 3 {x 2 ) and unknown lengths li, I2 such that 


u,2{h) = 7a = 0 , u^2{h) = 7 b = -2 cot 

I 3 {h) = I 3 a = 0 , iSik) = I 3 b = 2 cot ip. ( 5 ) 


Although the analytical solution can be found in the most general case with 
A 7^ 0 , we choose A = 0 for the shortness of presentation. In this special case 
/ (m, 2,/3) and g{( 3 ^ 2 ) are given by 


/ {u^2i ~ ^ 1 ^ sin^ (p + /d sin (p cos (p + l)‘ — 1 


+ (/9u^ 2 sin cos +/? cos^ — ^^2) + (/3sin(y9cos</9 — 1)^ — 1 

+2 [(—/dM,2 sin ip cos ip — 13 cos^ ip + ^^2) {( 3 u ^2 sin^ ip + (3 simp cos + l) 
n 12 'I 1/^9 sin^ ip 

+/?sinV(l-/^sinvjcosv?)] g (A2) = i^k ’ . 


Functional Q can be reduced to the functional depending only on [ 3 . Indeed, 
the variation of (|^ with respect to u{x2) leads to the equation and boundary 
conditions that imply 


du^2 


= 0 , 


so the stresses inside the boundary layer are zero. Since it was already 
established in [8] that the stresses in layers A and B are also zero, the whole 
specimen is stress-free. Solving the above equation with respect to u^2 we 
express it in terms of (3 


/d (/? sin 2(p — 2 cos 2 ip) 

(^cos2(p — 1) + 2/dsin2(p — 2’ 

Substituting this back into the functional Q we reduce the latter to the 
functional depending only on I 3 {x 2 ) and /i, I2 



\p{f 3 ) + g{f 3 , 2 )]dx 2 , 


( 6 ) 
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where 


p(/3) = f{u^2{P),/3) = P'^sm'^ip{/3smip -2cos(pf 

1 3/3^ - 8/33 sin 2(fi + 4/?3 sin 4^? - 4 (/3^ + 2) /3^ cos 2(fi 

1 8 [(3“^ (1 — cos 2lp) -\- (3“^ — 2(3 sin 2(^9 + 2]^ 

(/3^ — 4) (3"^ cos Alp + 12/3^ — 16/3 sin 2^9 + 16 1 
8 [(3‘^ (1 — cos 299 ) + (3“^ — 2(3 sin 2Lp + 2f‘ j 

Varying functional ([^ with respect to (3(x2), h, h and taking into account 
the constraints ([^ we obtain the Euler equation 

^g'(P.2)-p'W = o (7) 

( 1 X 2 

which is subjected to the boundary conditions 

(3{li) = Q, /l(/ 2 ) = 2cot(p, 

^2(^1) = 0 ) ^2(^2) = 0 . 

The last two boundary conditions mean that the dislocation density must 



Figure 3: Phase portrait and the separatrix connecting the states A and B. 

vanish at the boundaries X 2 = h and X 2 = h- Equation Q possesses the 
Erst integral 

g{(3p) - p{(3) = c. (9) 
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Using the boundary conditions (|^ as well as the identities p{I3a) = pif^s) = 
0, we hnd that c = 0. The phase portrait computed for b = 2.5 x 
ps = k = 10“® and ip = —7r/4 is shown in Fig. We see that 

the phase curve connecting the states A and B and satisfying the boundary 
conditions ([^ is the separatrix in the lower half of the phase plane (/9, (3^2) 
represented by the red curve with an arrow denoting the direction of change 
of (3 as X 2 changes from li to 12- It is also easy to show that the phase curve 
connecting B and A is the separatrix in the upper half of the phase plane, so 
the dislocations in the boundary layer between B and A have another sign 
than those connecting A and B. 



Figure 4: Plot of /3 versus the dimensionless coordinate ( 


y/2x2hps 

\/k\ sin ip\ 


The hrst integral ([^ with c = 0 enables one to hnd f3 implicitly in terms 
of X 2 according to 

■\/fc|sin(y9| dt 

^/^bps Jcotifi y/ p{t) 

The plot of versus the dimensionless coordinate ( = \/2x2bps/{'/k\ sin(p|) 
for the above chosen parameters is shown in Fig. It is seen that, for large 
ICI, the plastic slip remains in the very close neighborhood of the state A 
with /3 = 0 or the state B with (3 = 2 cot (p, while the strong change of (3 
leading to the transition from A to B occurs in a hnite interval of <(■ 

It turns out that if we compute the thickness of the boundary layer ac¬ 
cording to the formula 

y/k\ sin(p| dt 

^ — ^2 — — -/=- / — ! 1 

y2hps J2 cot ip yp{t) 
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Figure 5: Thickness of geometrically necessary boundary h (in nanometer) versus the 
misorientation angle 9 = 2ip + t: (in degrees). 

then this thickness becomes inhnite, so the result contradicts our assumption 
about the smallness of h as compared with sH, (1 — s)H, and with the 
sizes of the specimen. The resolution of this conflict should be found in the 
discreteness of crystals that may accommodate only a large but hnite natural 
number of dislocations. Since the dislocation density is p = |/9,2 sin (p|/6, the 
total number of dislocations in the boundary layer equals 



N = L 


As the smallest number of dislocation is 1, the smallest quantum of plastic 
slip leading to the recognizable change of N must be f3q = h/ (L| sin<p|). Now, 
if we take 6/L as a small positive number and compute the thickness of the 
boundary layer according to 



( 10 ) 


then the thickness becomes hnite (except for (p = 0). Since the pre-factor, 
the integrand, and both limits of integration depend on the orientation of slip 
system <p, the thickness is also a function of <p. Fig.j^shows the dependence of 
the thickness of geometrically necessary boundary (measured in nanometer) 
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on the misorientation angle 0 = 2^9 + tt (measured in degrees, see [S]) for 
the above chosen parameters and h/L = (the result turns out to be 

not strongly sensitive to the choice oi h/L). It is seen that the thickness 
changes from 50 to 6 nanometers for misorientation angles larger than 10°. 
For misorientation angles close to zero formula (10) is not well-dehned. In 
this case the dislocation model of small-angle tilt boundary considered in 
is more appropriate. 


YoCN/m) 



Figure 6: Energy density of of geometrically necessary boundary (in N/m) versus the 
misorientation angle 9 = 2ip + t: (in degrees). 

The minimum of the functional ([^ can easily be computed based on the 
analytical solution found above. Using the hrst integral we obtain 

Ib = 2 f p{l3)dx2 = / Vpil3)dl3. (11) 

Jh V ^OPs J2cot(/3 

Since this minimum value has the meaning of the energy (normalized by p) 
per unit area of the geometrically necessary boundary, the energy density 
of such boundary must be 7 g = ph- Similar to ( p)0| , the pre-factor, the 
integrand, and the lower limit of integration in ( [IT| ) depend on the orienta¬ 
tion of slip system p, so 7 g must also be a function of p. Fig. shows the 
dependence of 7 g (measured in N/m) on the misorientation angle 6 = 2p + 7r 
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(measured in degrees) for the above chosen parameters and /i = 26GPa (for 
aluminum). In contrast to h, the energy density of geometrically necessary 
boundary is well-dehned for all misorientation angles. For small misorien- 
tation angles up to 50° the energy density of the geometrically necessary 
boundary changes from zero to approximately 1 N/m that agrees quite well 
with the value 0.625 N/m given in [7j. 

In [8] we gave the estimation of the number of grains based on the follow¬ 
ing deliberations. The energy of the boundary layers near the side boundaries 
required to satisfy the side boundary conditions turns out to be of the order 
fiHDe, where £ is the thickness of one pair of layers A and B. Taking into 
account the energy of the geometrically necessary boundaries which is of the 
order 'jcDLH/e, we can estimate £ by minimizing these two contributions 
to the energy yielding 

e ~ \/jgLJJ1 ~ \/ hL. 

This relation exhibits clearly the size effect. 

4. Conclusion 

We have shown in this paper that the presence of the gradient of plastic 
slip in the energy of the dislocation network enables one to regularize the 
non-convex energy minimization in the class of smooth displacements and 
plastic slips. This leads to the formation of grains with geometrically nec¬ 
essary boundaries having a hnite thickness. Let us mention the similarity 
of the proposed continuum theory with the static version of Cahn-Hilliard 
theory of phase separation [3] and the recent phase-held approach [5]. The 
generalization of our model to capture the motion of geometrically necessary 
boundaries and the grain growth will be addressed elsewhere. 
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